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INTRODUCTION 


Paired  ionization  chambers  used  for  the  dosimetry  of 
mixed  gamma-neutron  radiation  consist  of  tissue- 
equivalent  (TE)  chambers  and  nonhydrogenous  cham¬ 
bers.  TE  chambers  respond  to  both  neutron  and  gamma 
radiation,  and  their  readings  are  closely  related  to  the 
total  dose.  Nonhydrogenous  chambers  are  relatively 
insensitive  to  neutron  radiation;  they  primarily  record 
the  gamma  dose.  This  report  presents  detailed  deriva¬ 
tions  of  the  characteristics  of  response  for  two  ioniza¬ 
tion  chambers  used  for  mixed  gamma-neutron  dosim¬ 
etry:  an  A-150  plastic  chamber  operated  with  TE  gas 
(TE/TE)  and  a  graphite  chamber  filled  with  carbon 
dioxide  (C/CO2).  (The  compositions  of  A-150  plastic- 
and  methane-based  TE  gas  are  given  in  reference  1.) 
These  spherical  50-cm^  ion  chambers  were  originally 
described  by  Lynn  (2)  and  Sayeg  (3),  and  are  dia¬ 
grammed  in  Figure  1.  They  served  as  the  backbone  of 
AFRRI  mixed-field  dosimetry  measurements  since  1965 
(e.g.,  references  4-6). 


-0.591t 
- — 1  265" - 


Figure  1.  AFRRI  spherical  50  cm^  ion  chambers 
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The  recent  availability  of  spectral  information  for  both 
neutron  and  gamma-ray  energy  (7,  8)  for  the  AFRRI 
TRIGA  reactor  has  prompted  the  present  review  of  the 
response  characteristics  of  the  50-cm^  chamber.  This 
information  enables  calculation  of  both  the  neutron  and 
gamma-ray  sensitivities  of  the  paired  chambers,  in 
accordance  with  currently  accepted  neutron  dosimetry 
protocols  (1,  9,  10).  Moreover,  the  physical  constants 
involved  in  calculating  neutron  dosimetry,  such  as 
kerma  factors  and  W-values,  are  known  now  with  much 
greater  certainty  (11-13).  Thus,  reevaluation  of  the 
response  characteristics  of  the  ion  chamber  is 
necessary,  to  serve  as  a  basis  for  future  work  and  to 
clarify  the  accuracy  of  previous  estimates  of  these 
parameters. 

Theory 

The  response  of  paired  ion  chambers  to  mixed  neutron- 
gamma  radiation  is  expressed  as  follows  (1): 

Rt  =  k^  Dj<r  +  h'p  Dq  (Eq  1) 

Ru  =  ku  dN  +  hU  dG  (E<1  2) 

where  Djj  and  Dq  represent  the  respective  neutron  and 
gamma  components  of  the  total  dose  deposited  in 
tissue.  Rt  and  Ry  represent  the  readings  of  TE/TE  and 
C/CO2  ion  chambers,  respectively,  divided  by  the 
cobalt-60  calibration  factor  (coulomb  per  rad)  for  each 
chamber.  The  k's  and  h's,  known  as  the  paired  chamber 
constants,  represent  the  sensitivity  of  each  chamber 
(subscript  T  for  TE/TE,  U  for  C/CO2)  to  neutron  and 
gamma  radiation  (k  for  neutron,  h  for  gamma),  relative 
to  its  sensitivity  to  cobalt-60  radiation. 
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Values  of  the  paired  chamber  constants  depend  on  the 
energy  spectra  of  the  neutron  and  gamma  radiation, 
according  to  the  following: 


WC  ,  .  (hen/p)tApen/p)nn  q 

k-r  1 1  =  ' — -  •  (Smg)p 

'  WN  L  (Kt/Kg)N 


(Eq  3) 


Wq  (Smg>Q  (Pen/plt^Men/p)111  Q 

h  j  (j  =  — “  * - 

WG  (Smg)G  (Pen/p)t/(lJen/p)m  G 


(Eq  4) 


where 


C  = 
N  = 

G  = 

m  = 

g  = 
t 

W  = 
(Smg)  = 


( y  en/  P )  - 
K  = 


Subscript  for  cobalt-60  calibration  field 
Subscript  for  neutron  spectrum  to  be 
measured 

Subscript  for  gamma  spectrum  to  be 
measured 

Subscript  for  ion  chamber  wall  material 
Subscript  for  ion  chamber  gas 
Subscript  for  tissue,  i.e.,  ICRU  muscle  (1) 
Average  energy  to  produce  an  ion  pair  in 
TE  or  CO2  gas 

Ratio  of  electron  mass-stopping  powers  of 
ion-chamber  wall  material  to  that  of 
ion-chamber  gas  material 
Mass  energy  transfer  coefficient  for 
gamma  radiation 

Average  kerma  factor  of  neutron 
spectrum 
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That  the  above  equations  apply  to  the  AFRRI  50-cm^ 
paired  ion  chambers  is  based  on  two  assumptions: 
(a)  Ionizations  in  chamber  gas  caused  by  neutrons  are 
entirely  due  to  neutron  interactions  that  occur  in  the 
gas,  not  the  wall.  This  is  because  the  ion  chambers  are 
relatively  large,  compared  to  the  short  range  of  recoil 
nuclei  generated  in  neutron  interactions  in  chamber 
walls.  This  is  especially  true  of  the  low-energy 
neutrons  encountered  in  fission  neutron  spectra.  For 
C/CO2  chambers,  Markewicz  and  Pszona  (14)  showed 
that  even  a  1-cm^  ion  chamber  satisfied  this  condition 
for  neutrons  with  energies  up  to  15  MeV.  For  TE/TE 
chambers,  recoil  protons  from  chamber  wall  material 
have  longer  ranges  than  do  recoil  C  nuclei  in  the  C/CO2 
chamber,  but  the  kerma  factors  for  A-150  plastic  and 
TE  gas  are  almost  equal  to  each  other  (1),  and  the 
W-value  for  wall  recoils  is  nearly  equal  to  the  W-value 
for  gas  recoils  (12).  (b)  The  above  equation  for  hj  and 
hy  uses  Fano's  Theorem,  based  on  the  nearly  matching- 
elemental  composition  of  wall  and  gas  in  the  TE/TE 
chamber,  and  on  the  similarity  in  atomic  numbers  of 
carbon  and  oxygen,  together  with  the  almost  unity 
value  of  (Smg)  (see  below)  for  the  C/CO2  chamber. 

METHODS 

Cobalt-60  Calibration  Field 


Dosimetric  parameters  for  the  cobalt-60  calibration 
field  are  given  in  Table  1,  with  appropriate  references. 
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Table  1 .  Dosimetric  Parameters  for  Cobalt-60  Calibration  Field 


Parameter 

Value 

Reference 

(y  en/P  M^en/P  )m 

TE/TE 

1.001 

Broerse  et  al.  (1981) 

c/co2 

1.126 

ORNL  (1975) 

(Smg)c 

TE/TE 

1.000 

Broerse  et  al.  (1981) 

C/C02 

1.009 

AAPM  (1980) 

Wc 

TE/TE 

29.3  JCT1 

Broerse  et  al.  (1981) 

c/co2 

33.0  JC_1 

1CRU  (1979) 

Gamma  Fields  of  the  Reactor 


Gamma  fields  of  the  reactor  include  photons  distributed 
in  energy  from  10  keV  to  10  MeV  (8).  The  average 
gamma  energies  of  these  spectra  range  from  0.67  to  1.8 
MeV  (15).  For  such  spectra,  the  average  energy  to 
produce  an  ion  pair,  W,  can  be  taken  to  be  the  same  as 
for  cobalt-60  gamma  rays  (16). 

For  the  TE/TE  ion  chamber,  the  tissue  equivalence  of 
the  material  of  the  chamber  wall  and  the  close  match 
in  elemental  composition  of  the  wall  and  gas  allow  use 
of  the  same  stopping  power  ratio  and  mass  energy 
transfer  coefficient  ratio  for  reactor  gamma  fields  as 
for  cobalt-60  gamma  rays.  Thus,  for  the  TE/TE  ion 
chamber,  the  paired  chamber  constant  h-p  can  be  taken 
to  be  unity  for  all  reactor  gamma  fields,  as  is  done  for 
higher  energy  neutron  beams  (9,  10). 
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For  the  C/CO2  ion  chamber,  the  energy  dependence  of 
both  (Smg)  and  y  en/  p  must  be  accounted  for  to 
properly  evaluate  hyj  for  energy  spectra  of  reactor 
gamma  rays.  This  was  completed  in  an  earlier  report 
(15)  for  the  mass  energy  transfer  coefficients;  results 
are  duplicated  in  Table  2.  Because  the  gas-to-wall 
ratio  (Smg)  for  mass-stopping  power  is  a  slowly  varying 
function  of  electron  (or  gamma-ray)  energy,  an  approxi¬ 
mate  method  was  used  to  arrive  at  suitable  values  for 
reactor  gamma  fields.  For  each  reactor  configuration, 
the  average  gamma-ray  energy  was  used  to  derive  an 
average  energy  for  Compton  recoil  electrons  (17).  This 
average  energy  was  referred  to  tables  of  electron 
collision  mass-stopping  power  (18),  to  derive  the  ratios 
of  stopping  power  shown  in  Table  2.  These  ratios  varied 
between  1.002  and  1.022,  compared  to  1.009  for  cobalt- 
60  gamma  rays.  The  lowest  ratios  were  for  the 
configurations  with  the  largest  populations  of  low- 
energy  scattered  photons,  i.e.,  the  configurations  near 
the  back  walls  of  the  exposure  rooms. 

Reactor  Neutron  Fields 

The  ratios  of  spectrum-averaged  kerma  factors  for 
tissue  to  those  for  both  TE  gas  and  CO2  were  derived  in 
an  earlier  report  (15),  and  are  listed  in  Table  3  for 
reference.  Neutron  W-values  for  TE  gas  and  CO2  have 
been  the  subject  of  several  recent  reports  (12,  13) 
demonstrating  a  complicated  dependence  on  energy  for 
both  gases.  Published  data  on  neutron  W-values  is 
incomplete,  particularly  for  low-energy  neutrons  (below 
0.5  MeV),  and  the  remaining  uncertainties  in  W  have 
been  recognized  as  one  of  the  major  sources  of 
uncertainty  in  dosimetry  of  neutron  radiation  (1).  The 
W-values  used  in  the  present  report  are  given  in  Table  4 
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for  the  energy  groupings  in  which  reactor  neutron 
spectra  were  tabulated.  Detailed  discussion  of  these 
values  is  presented  in  the  Appendix.  Averaging  of  the 
neutron-W  values  over  the  reactor  neutron  spectra  was 
arrived  at  by  using  computer  program  SPT  (15),  as 
described  in  the  Appendix.  Resulting  spectrum- 
averaged  W-values  for  all  reactor  configurations  are 
given  in  Table  3. 


Table  2.  Dosimetric  Parameters  for  Graphite-C02  Ion  Chamber  in  Gamma 
Ray  Fields  of  Reactor 


Room 

Distance 
to  Core 
(cm) 

Configuration* 

Average 
Energy  (MeV) 

(Smg)Q 

(^en/PV 

(Uen/P^mQ^ 

1 

50 

Unshielded 

1.16 

1.012 

1.128 

100 

Unshielded 

1.15 

1.012 

1.130 

200 

Unshielded 

1.05 

1.010 

1.132 

300 

Unshielded 

0.91 

1.007 

1.134 

400 

Unshielded 

0.78 

1.004 

1.139 

500 

Unshielded 

0.67 

1.002 

1.149 

2 

50 

Unshielded 

1.18 

1.012 

1.129 

100 

Unshielded 

1.13 

1.011 

1.129 

200 

Unshielded 

0.95 

1.008 

1.135 

300 

Unshielded 

0.76 

1.004 

1.142 

2 

30 

Pneumatic  tubes 

0.83 

1.005 

1.138 

1 

100 

5"  H20 

1.04 

1.009 

1.131 

100 

12"  H20 

0.95 

1.009 

1.137 

100 

2"  Pb 

1.65 

1.008 

1.128 

100 

2"  Pb  and  exercise 
wheel 

1.36 

1.022 

1.131 

100 

6"  Pb 

0.92 

1.016 

1.142 

100 

6"  Pb  and  cave 

1.80 

1.007 

1.132 

100 

6"  Pb  and  cylinder 
phantom 

1.13 

1.011 

1.133 

*For  complete  description  of  reactor  configuration,  see  references  7  and  8. 
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Table  3.  Dosimetric  Parameters  for  Neutron  Fields  of  Reactor 


Room 

Distance 
to  Core 
(cm) 

Average  Neutron 

Configuration  Energy  (MeV) 

WN  (eV) 
TEgas  C02 

Kt/Kg 

TEgas  co2 

1 

50 

Unshielded 

1.46 

31.7 

44.3 

0.980 

9.57 

100 

Unshielded 

1.49 

31.7 

44.2. 

0.979 

9.50 

200 

Unshielded 

1.32 

31.6 

44.4 

0.983 

9.62 

300 

Unshielded 

1.16 

31.7 

44.5 

0.981 

9.89 

400 

Unshielded 

1.11 

31.8 

44.6 

0.986 

9.78 

500 

Unshielded 

0.89 

31.8 

45.2 

0.984 

9.84 

2 

50 

Unshielded 

1.35 

31.7 

44.1 

0.979 

9.59 

100 

Unshielded 

1.32 

31.7 

44.4 

0.986 

9.66 

200 

Unshielded 

1.11 

31.7 

44.2 

0.984 

9.68 

300 

Unshielded 

0.93 

31.7 

44.8 

0.981 

9.81 

2 

30 

Pneumatic  tubes 

0.88 

31.8 

44.9 

0.982 

10.21 

1 

100 

5"  H20 

1.99 

31.5 

42.9 

0.979 

9.26 

100 

12"  H20 

3.71 

31.2 

40.7 

0.983 

8.49 

100 

2"  Pb 

1.58 

31.7 

45.0 

0.986 

9.67 

100 

2"  Pb  and  exercise 
wheel 

1.05 

31.8 

45.7 

0.984 

9.92 

100 

6"  Pb 

0.96 

32.0 

47.1 

0.983 

10.18 

100 

6"  Pb  and  cave 

0.68 

32.0 

49.0 

0.986 

10.38 

100 

6"  Pb  and  cylinder 
phantom 

0.55 

32.3 

47.7 

0.986 

11.00 

Table  4. 

Group  W  Values  of  the  Neutron  Spectru 

Energy 

Group  Upper 

W(eV) 

Group 

Energy  (MeV) 

TE  Gas 

co2 

1 

19.6 

31.0 

35.5 

2 

16.9 

31.0 

35.9 

3 

14.9 

31.0 

36.1 

4 

14.2 

31.0 

36.1 

5 

13.8 

31.0 

36.4 

6 

12.8 

31.0 

36.7 

7 

12.2 

31.0 

37.0 

8 

11.1 

31.1 

37.5 

9 

10.0 

31.1 

37.9 

10 

9.0 

31.1 

38.4 

11 

8.2 

31.2 

39.0 

12 

7.4 

31.0 

38.7 

13 

6.4 

31.1 

40.1 

14 

5.0 

31.1 

39.3 

15 

4.7 

31.2 

40.2 

16 

4.1 

31.5 

42.4 

17 

3.0 

31.3 

41.7 

18 

2.4 

31.0 

41.7 

19 

2.3 

31.3 

43.8 

20 

1.8 

31.5 

47.1 

21 

1.1 

32.9 

49.0 

22 

0.55 

32.4 

61.3 

23 

0.16 

33.4 

75.2 

24 

0.11 

33.6 

84.1 

25 

0.052 

34.6 

98.2 

26 

0.025 

36.0 

107.7 

27 

0.022 

38.3 

28 

0.010 

51.0 

29 

0.0034 

146.2 

30 

0.0011 
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RESULTS 


Before  presenting  the  derived  paired-chamber  constants 
for  the  various  reactor  configurations,  each  ion 
chamber's  energy  response  curves  for  both  gamma  and 
neutron  radiation  should  be  examined.  The  neutron 
results  (shown  in  Figure  2)  were  obtained  using  the 
neutron-W  values  of  Table  4  and  the  neutron-kerma 
factors  (11)  for  individual  spectrum-energy  groups. 
Figure  2  also  shows  the  response  curve  calculated  by 
Markewicz  and  Pszona  (14),  using  Monte  Carlo  tech¬ 
niques  for  a  large  C/CO2  chamber  for  neutrons  above  1 
MeV.  This  curve  compares  closely  with  the  results 
derived  here.  The  range  of  reported  C/CO2  kjj  values 
summarized  by  Mijnheer  (19)  for  15-MeV  neutrons  also 
compares  favorably  with  the  present  results.  Measured 
TE/TE  (1  cm3)  and  C/CO2  (2  cm^)  ion-chamber 
responses  for  neutrons  over  1  MeV  reported  by 
Waterman  et  al.  (20)  are  also  shown  in  Figure  2.  For 
the  TF/TF  chamber,  the  present  results  agree  well  with 
the  measured  responses  for  neutrons  below  10  MeV. 
The  lower  C/CO2  ky  values  measured  for  neutrons 
above  5  MeV  is  surprising,  because  for  smaller  ion 
chambers  (such  as  those  used  by  Waterman  et  al.),  one 
would  expect  ky  values  to  be  higher  than  for  large 
chambers  (14).  Possible  explanations  for  the  low  k^p  and 
ky  values  measured  at  high-neutron  energies  are 
discussed  by  Waterman  et  al.  (20). 

At  the  time  of  this  writing,  for  neutron  energies  below 
1  MeV,  there  exists  no  published  calculation  or 
measurements  of  kj  or  ky  to  compare  with  the  present 
results.  This  is  due  to  the  emphasis  on  high-energy 
neutron  beams  in  the  recent  literature,  particularly  for 
clinical  radiotherapy  applications. 
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Neutron  Energy  (MeV) 

Figure  2a.  Neutron  sensitivity  of  TE  ion  chamber 


The  calculated  photon  energy  response  of  the  C/CO2 
ion  chamber  is  shown  in  Figure  3.  (Recall  that  the 
photon  energy  response  of  the  TE/TE  ion  chamber  was 
taken  to  be  unity,  based  on  the  tissue  equivalence  of 
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A-150  plastic  and  methane-based  TE  gas.)  The 
decrease  in  hy  for  the  C/CO2  ion  chamber  at  low 
photon  energies  occurs  mainly  because  the  graphite 
atomic  number  is  lower  than  that  of  tissue,  thus  causing 
fewer  photoelectric  interactions.  At  high  photon  ener¬ 
gies  (over  2  MeV),  the  increased  values  of  hy  are  due 
mainly  to  the  energy  dependence  of  (Smg).  Recall  that 
(Smg)  was  obtained  by  an  approximation  technique  (see 
Methods),  but  the  relatively  low  flux  of  reactor  gamma 
rays  above  2.2  MeV  minimizes  errors  in  the  photon 
response  due  to  (Smg)  uncertainties. 


Figure  3.  Gamma  sensitivity  of  C/COj  ion  chamber 


Paired  ionization  chamber  constants  evaluated  for 
TRIGA  reactor  gamma-neutron  spectra  are  given  in 
Table  5.  The  precision  of  displayed  results  is  probably 
greater  than  warranted  by  the  accuracy  of  the  calcula¬ 
tions,  but  it  is  useful  for  demonstrating  the  spectral 
dependence  of  the  constants.  These  results  show  a 
clear  separation  in  kj  and  ky  values  between  the 
water-shielded  reactor  configurations  and  all  other  con¬ 
figurations.  This  occurs  because  of  the  higher  average 
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energies  of  water-shielded  neutron  spectra.  For  hy,  a 
single  value  was  found  to  apply  within  ±  1%  for  all 
spectra. 


able  5. 

Paired  Ion 

Distance 
to  Core 
(cm) 

Chamber  Constants 

for  AFRRI 

TE/TE 

Reactor  Configurat 

c/co2 

Room 

Configuration 

kT 

ku 

hu 

1 

50 

Unshielded 

0.944 

0.0862 

0.995 

100 

Unshielded 

0.945 

0.0870 

0.993 

200 

Unshielded 

0.944 

0.0855 

0.994 

300 

Unshielded 

0.943 

0.0830 

0.995 

400 

Unshielded 

0.935 

0.0837 

0.994 

500 

Unshielded 

0.937 

0.0821 

0.987 

2 

50 

Unshielded 

0.945 

0.0864 

0.994 

100 

Unshielded 

0.938 

0.0852 

0.995 

200 

Unshielded 

0.940 

0.0854 

0.993 

300 

Unshielded 

0.943 

0.0831 

0.991 

2 

30 

Pneumatic  tubes 

0.939 

0.0797 

0.993 

1 

100 

5"  H20 

0.951 

0.0919 

0.996 

100 

12"  H20 

0.956 

0.1057 

0.990 

100 

2"  Pb 

0.938 

0.0839 

0.999 

100 

2"  Pb  &  exercise 
wheel 

0.937 

0.0806 

0.983 

100 

6"  Pb 

0.933 

0.0762 

0.979 

100 

6"  Pb  &  cave 

0.930 

0.0718 

0.997 

100 

6"  Pb  <5c  cyl 
phantom 

0.921 

0.0696 

0.992 

*hT  =  1.000  for  all  configurations. 


DISCUSSION 

This  work  has  incorporated  recently  reported  AFRRI 
TRIG  A  reactor  gamma-neutron  energy  spectra  together 
with  up-to-date  estimates  of  pertinent  kerma  factors 
and  W-values,  to  evaluate  the  energy  response  charac¬ 
teristics  of  AFRRI  50-cm^  paired  ionization  chambers. 
This  derivation  of  the  paired-chamber  constants  shall 
serve  as  the  basis  for  future  dosimetric  measurements 
with  these  ion  chambers;  it  can  also  elucidate  the 
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validity  of  estimates  of  these  constants  used  in  the 
past.  Verification  of  the  derived  constants  by  experi¬ 
mentation  has  not  been  considered  in  the  present  work, 
but  it  will  be  the  subject  of  future  research. 
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APPENDIX. 


VALUES  OF  NEUTRON- W  FOR  METHANE-BASED 
TE  GAS  AND  CO2 
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W  is  defined  as  the  quotient  of  E  divided  by  N,  where  N 
is  the  mean  number  of  ion  pairs  formed  when  the  initial 
kinetic  energy,  E,  of  a  charged  particle  is  completely 
dissipated  in  the  gas.  For  neutron  radiation,  it  is 
necessary  to  consider  the  energy  deposited  by  several 
different  types  of  charged  particles  interacting  in  the 
gas.  Following  Goodman  and  Coyne  (12),  the  W-value 
for  neutrons  of  energy,  En,  is  defined  as  follows: 


WN(EN) 


f  Ki<ENl 
?  Nj(EN) 


(Eq  A-J) 


Where  Nj(En)  is  the  number  of  ions  produced  by  second¬ 
ary  charged  particles  of  type  j,  which  deposit  kerma 
Kj(En).  Nj(En)  and  Kj(En)  are  calculated  as  follows: 


-max  j 


N  j  ( E  N ) 


j  '  nj(E,EN)EdE 


Wj(E) 


(Eq  A-2) 


Emax  j 

Kj(EN)  =  j  nj ( E , E N ) E d E 
o 


(Eq  A-3) 


where  nj(E,En)  represents  the  initial  differential  energy 
spectrum  of  charged  particles  of  type  j  produced  by 
unit  neutron  fluence  of  energy  En,  and  W7j(E)  is  the 
W-value  for  type  j  particles  with  initial  energy  E. 


For  neutrons  with  energies  distributed  over  a  spectrum, 
the  above  theory  must  be  extended  to  define  a 
spectrum-averaged  neutron  W-value,  denoted  W*n,  as 
follows: 


WN  =  Ed/Ns 


(Eq  A-4) 
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In  the  above  equation,  E<j  is  the  energy  deposited,  and 
Ns  is  the  total  ions  produced  per  unit  fluence.  E<j  and 
Ns  are  determined  as  follows: 


00 


Ed  =  )  K(EN)0(EN)dEN  =  2^  Kg0gAEg  (Eq  A-5) 


o 


CO 


where  K(En)  represents  the  gas-kerma  factor  for 
neutrons  of  energy  En,  and  Kg  represents  the  gas-kerma 
factor  for  neutrons  that  are  in  group  g.  Then  <f>  (E)  and 
(pg  represent  the  neutron  energy  spectrum,  with  the 
sums  extending  over  the  37  energy  groups  in  the  neu¬ 
tron  spectral  format  of  reference  21.  To  complete  the 
sums  in  the  above  equations,  it  is  necessary  to  use 
average  Wn  values  for  each  energy  group  (denoted  by 
Wg).  These  are  derived  by  averaging  Wn  between  the 
upper  Ey  and  lower  El  energies  of  each  group,  as 
follows: 


(Eq  A- 7) 


W, 


9 


To  obtain  the  right  side  of  equation  A-7,  the  neutron 
kerma  factor  K(En)  has  been  assumed  to  be  proportional 
to  En  between  neutron  energies  El  and  Ey,  and  <j>(En) 
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within  each  energy  group  has  been  assumed  to  have  a 
dependence  on  1/En. 

TE  GAS 

The  neutron  W-values  for  TE  gas  of  Goodman  and 
Coyne  (12)  were  used  as  the  basis  for  the  W-values  of 
the  spectrum  groups  shown  in  Table  4  for  the  energy 
range  19.6  MeV  through  110  keV  (energy  spectrum 
groups  1  through  23).  In  the  following,  W  is  estimated 
for  TE  gas  for  neutrons  with  energies  between  1  and 
110  keV  (groups  24  through  29).  These  estimates  were 
based  on  the  following  assumptions: 

a.  The  neutron  kerma  and  the  ionizations 
produced  by  neutrons  in  TE  gas  were  con¬ 
sidered  to  be  due  to  recoil  protons  for  the 
purpose  of  evaluating  equations  A-l,  A-2,  A-3, 
and  A-7.  This  assumption  was  based  on  the 
fact  that  9496-96%  of  the  total  kerma  is,  in 
fact,  due  to  recoil  protons  in  this  range  of 
neutron  energy,  and  this  percentage  varies 
only  slowly  with  neutron  energy.  The 
W-values  calculated  on  the  basis  of  this 
assumption  were,  in  the  end,  corrected  for  the 
fraction  of  total  kerma  due  to  C,  N,  and  O. 

b.  Goodman  and  Coyne's  (12)  empirical  expres¬ 
sion  for  W  for  protons  in  TE  gas,  which  is  valid 
for  proton  energies  above  4  keV,  was  extrap¬ 
olated  to  lower  proton  energies.  This  expres¬ 
sion  is  Wp  =  29.28  (in  E)-^  +  29.99,  where 
proton  energy  E  is  in  units  of  keV  and  Wp  is  in 
units  of  eV.  This  assumption  was  needed 
because  no  experimental  data  existed  on  W  for 
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protons  (or  neutrons)  with  energies  below  4 
keV.  Extrapolation  of  this  expression  to  lower 
energies  had  no  physical  justification,  but  no 
viable  alternative  was  available. 

c.  Recoil  protons  with  energies  below  1  keV  were 
assumed  to  produce  no  ionizations.  Thus,  the 
lower  integration  limit  of  equation  A-2  for  Nj 
is  set  at  1  keV. 

Based  on  the  above  assumptions,  the  initial  differential 
energy  spectrum  np(E,En)  for  recoil  protons  is  given  by 
the  spectrum  that  describes  elastic  scattering,  namely 
np(E,En)  =  <t>  p(En)/En,  where  c}>p(En)  represents  the 
number  of  recoil  protons  arising  from  elastic  scattering 
with  neutrons  of  energy  En.  Equations  A-2  and  A-3 
were  evaluated  as  follows,  with  the  subscript  p  denoting 
recoil  protons: 


Kp(EN)  =  %EN*p(EN) 


(Eq  A-8) 


(Eq  A- 9) 


The  above  integral  for  Np(Ep)  was  evaluated  numer¬ 
ically,  and  then  approximated  by  the  following 
expression  with  an  error  of  3%  at  3  keV  and  less  than 
1%  at  >  4  keV: 


(Eq  A-10) 
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Using  the  results  of  equations  A-8  and  A-10  in  equation 
A-l  gives  the  following  for  neutron  W-values  in  TE  gas: 


Wn(En)  =  30.86 


(Eq  A-ll) 


Recall  that  the  above  result  does  not  include  a  correc¬ 
tion  for  the  fraction  of  the  neutron  kerma  due  to  C,  N, 
and  O  in  TE  gas.  Equation  A-ll  was  integrated 
following  equation  A-7,  to  arrive  at  neutron  Wg-values 
for  energy  spectrum  groups  bounded  by  the  lower  and 
upper  energies  El  and  Ey,  respectively. 


(Eq  A-l 2) 


In  the  above  equation,  the  factor  f  represents  the 
fraction  of  the  total  neutron  kerma  in  TE  gas  that  is 
due  to  recoil  protons.  (From  reference  11,  f  is  equal  to 
0.943  for  energy  group  23,  0.950  for  24,  0.956  for  25, 
0.960  for  26  and  27,  0.962  for  28,  and  0.960  for  29.) 
Equation  A-l 2  yields  the  group  Wg-values  shown  in 
Table  4  for  spectrum  groups  24  through  29. 

For  spectrum  group  23  (110-160  keV),  equation  A-12 
yields  a  Wg-value  of  33.37  eV,  which  agrees  within  2.4% 
with  the  value  (32.60)  estimated  for  group  23  from  the 
results  obtained  by  Goodman  and  Coyne  (12).  That  the 
Wg-value  derived  here  is  higher  than  the  Wg-value  of 
Goodman  and  Coyne  is  probably  due  to  the  fact  that  the 
present  calculation  ignores  ions  produced  by  C,  N,  and 
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O  recoils.  For  neutron  energies  below  100  keV,  the 
percentage  of  ions  produced  by  C,  N,  and  O  is  even 
lower  than  in  spectrum  group  23,  so  errors  due  to  this 
omission  should  be  below  2%  for  spectrum  groups  24 
through  29. 

Carbon  Dioxide 


Neutron  W-values  from  three  sources  were  used  to 
arrive  at  the  spectrum  group  Wg-values  shown  in  Table 
4:  (a)  Rubach  and  Bishell's  (13)  19  data  points  between 
0.45  and  14.0  MeV  for  infinitely  large  ion  chambers 
(13).  (Our  graphite/C02  ion  chamber  is  50  cc  volume.) 

(b)  Dennis'  data  at  14.0,  14.7,  and  15.5  MeV  normalized 
to  Rubach  and  Bishell's  W revalue  at  14.0  MeV  (22). 

(c)  The  estimated  W-values  for  low-energy  neutrons 
derived  below.  These  Wj^-values  are  plotted  in  Figure 
4.  Following  equation  A-7,  interpolation  between 
1/W-values  at  discrete  energies  was  used  to  evaluate 
the  spectrum  group  Wg-values  listed  in  Table  4. 

W-values  for  low-energy  neutrons  (below  0.5  MeV)  in 
CO2  were  estimated,  based  on  the  following 
assumptions: 

a.  Neutron  interactions  are  dominated  by  elastic 
scattering  with  C  and  O.  Initial  differential 
energy  spectra  for  C  and  O  recoils  are  approx¬ 
imated  by  those  for  isotropic  scattering  in  the 
center  of  mass  system,  namely:  <f>c(En)/Ec(En) 
and  (j>0(En)/E0(En),  respectively,  where  Ec  = 
0.284  En  and  E0  =  0.221  En  are  the  maximum 
recoil  energies. 

b.  Neutron  W-values  for  recoil  C  and  O  ions  in 
CO2  are  taken  from  the  empirical  expressions 
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of  Rubach  and  Bischell  (13).  These  are  given 
below  for  recoil  energies,  E,  expressed  in  MeV. 


Wc  = 


(  6.33  -  15.91  In  E, 
I  38.26  -  3.26  In  E, 


E  <  0.080  MeV 
0.080  <  E  <  6  MeV 


(  4.57  -  19.5  In  E,  E  <  0.1  MeV 

)  45.59  —  2.99  In  E,  0.1  ^E  <3.5  MeV 


(Eq  A-13) 


Based  on  the  assumption  of  isotropic  elastic  scattering, 
the  kerma  integrals  (equation  A-3)  are  evaluated  as 
follows: 


Ko  “  %E o0o(eN) 
Kc  =  1/2Ec0c(EN) 


(Eq  A-14) 


Using  Rubach  and  Bischell's  (13)  empirical  expressions 
for  W,  equation  2  can  be  evaluated  for  the  quantities 
Nj(En).  Consider  first  the  C  recoils. 


0c  (E  i\ 


Nc(En)  = 


/ 


EdE 


-c  o 


6.33  -  15.91  In  E 


E  <  0.080  MeV 


(Eq  A-15) 


f>  c(EN) 


0.08 

/ 


EdE 


Ec 

/ 


EdE 


o  6.33  -  15.91  In  E  0os  38.26  -  3.26  In  E 


0.080  <  E  <  6  MeV 
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Both  the  C  and  the  O  expressions  for  Nj(En)  involve 
integrals  of  the  following  form,  which  by  a  transforma¬ 
tion  of  variables  are  reduced  to  the  Exponential 
Integral  Ei(t),  as  below: 


EdE 

a  +  P  In  E 


exp(-2a/j3) 

~P 


/ 


t  (a) 


e't  dt 
t 


-2 

where  t  =  —  (a  4*  j3  In  E) 

P 


(Eq  A-16) 


exp(— 2a//3) 

-P 


a))  -  Ei  (t(b)) 


] 


The  integrals  in  the  expressions  for  Nc(En)  of  equation 
A-15  were  evaluated  for  selected  neutron  energies, 
using  tables  of  the  Exponential  Integral  (23).  Table  6 
lists  the  results  of  these  evaluations  expressed  in  terms 
of  the  quantity  =  (EC2/2I),  which  is  effectively  the 
W-value  for  a  C-recoil  spectrum,  with  maximum  energy 
Ec.  Similarly,  the  quantity  N0(En)  was  evaluated  for 
oxygen  recoils,  and  terms  =  (Eg/21)  are  listed  in 

Table  6.  Using  these  quantities,  Nc(En)  and  N0(En)  can 
be  expressed  as  follows: 


NC(EN)=%EC0C(EN)£2C-1 
Nq  (E[\j)  =1/2Eo0o(E|\j)S2o  1 


(Eq  A-17) 


Combining  equations  A-14  and  A-17  gives  the  Wn- 
values  for  CO2  as 


W 


N 


1  +  f' 

fV1  +  f'  ^c'1 


(Eq  A-18) 


26 


where  the  quantity  f'  is  defined  as: 

f-  =  E.c0c!!n1  (Eq  A-19) 

Thus  f '  represents  the  fraction  of  the  neutron  kerma  in 
CO2  due  to  Orecoils  divided  by  the  fraction  due  to 
O-recoils.  This  ratio  is  defined  by  the  C  and  O  kerma 
factors  and  the  relative  abundance  of  C  and  O  in  CO2. 
Values  of  f '  calculated  from  data  in  reference  11  are 
given  in  Table  6.  The  data  in  Table  6,  together  with 
equation  A-18,  formed  the  basis  for  the  low-energy  Wn- 
values  (below  0.5  MeV),  plotted  in  Figure  4. 


Table  6.  Quantities  Used  in  Calculation  of  W  for  Low-Energy 
Neutrons  in  CO2 


En 

(keV) 

Ec 

^  c 

Eo 

^  0 

f' 

22 

6.2 

96.49 

4.9 

121.87 

.79 

25 

7.1 

93.70 

5.5 

118.69 

.79 

52 

14.8 

80.77 

11.5 

101.53 

.76 

110 

31.2 

68.66 

24.3 

86.25 

.72 

160 

45.4 

62.65 

35.4 

78.73 

.48 

445 

126 

48.31 

98.3 

58.38 

.27 

550 

156 

47.22 

122 

54.83 

.80 

*Eq  and  Eq  are  maximum  energies  of  C  and  O  recoil  nuclei,  respectively,  arising 
from  elastic  scattering  with  neutrons  of  energy  En.  and  are  respective 
W-values  for  recoil  C  and  O  nuclei,  with  energies  covering  spectrum  from  zero  to 
Eq  and  Eq,  respectively.  F  is  fraction  of  neutron  kerma  in  CO2  due  to  C-recoils, 
divided  by  fraction  due  to  Orecoils. 
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For  445-keV  neutrons,  the  above  derivation  gave  a 
Wn-value  of  55.90  eV,  which  differed  by  3.3%  from  the 
Wn-value  published  by  Rubach  and  Bischell  (13).  At  this 
energy,  the  neutron  cross  section  for  oxygen  shows  a 
strong  resonance  peak,  and  the  assumption  of  isotropic 
elastic  scattering  is  at  its  weakest. 
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